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Abstract 

A new method named rational expansion method of exponent function 
is presented to find exact traveling wave solutions of differential- difference 
equations. This method generalizes the so-called tanh-method and other 
similar methods. Some examples are dealt with and their abundant ex- 
act solutions which include solitary solutions and periodic solutions are 
obtained. Among them, many solutions are new. 
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1 Introduction 

In some physical fields such as solid physics and biophysics, difFerential-diflerence 
equations (DDEs) are used to model some physical phenomena ([Illl]) such as 
particle vibrations in lattices, currents in electrical networks, pulses in biologi- 
cal chains, etc. One also use DDEs to simulate numerically soliton dynamics in 
high energy physics and other problems where they arise as approximations of 
continuum models. DDEs are semi-discretized with some (or all) of their spe- 
cial variables discretized while time is usually kept continuous. Many aspects of 
DDEs such as integrability criteria, the computation of densities and symmetries 
and so on have been investigated([31 01 [51 [5]). Some works have been done to 
investigate the exact solutions of DDEs. Qian ct al(T) have extended multilin- 
ear variable separation approach to a special DDE. In particular, some symbolic 
computation methods have been applied to this problem. For example, recently 
some expansive methods such as tanh-method([5]) and its generalizations ([Hj) 



and elliptic function expansive method ([TO]) have been proposed to solve the 
traveling wave solutions of DDEs. Using symbolic computation software such 
as Maple or mathematical some exact solutions of some DDEs have been ob- 
tained. In present paper, we propose a new method named rational expansion 
method of exponent function for solving traveling wave solutions of DDEs. Our 
method unify tanh-method and its some generalizations. For illustration, we ap- 
ply our proposed method to some examples such as Langniuir lattice, discrete 
mKdV lattice equation and Hybrid lattice equation, and derive their abundant 
exact solutions. Among them, many solutions are new. Those solutions include 
abundant solitary solutions and periodic solutions, thus they will be helpful for 
physical investigation further. 

2 Rational expansion method of exponent func- 
tion 

Let us consider rather general lattice equation formed as: 

P{u'^{x),u'^{x), - ■ ■ ,ull"^^\x),Un+kAx)r ■ ■ ,Un+k^^{x)) =0, (1) 

where P is a polynomial of its entries, u, x and n all represent multi-components, 
and u^^^ denotes the collection of mixed derivative terms of order r. We take 
traveling wave transformation 



e = d,n, + J2 CjXj + (5, (2) 



where di,Cj and 8 are constants to be determined. Under the transformation, 
the equation (1) becomes 

where the prime is derivative with aspect to ^. A crucial step of our method is 
to assume that the solutions of the equation (3) has the following form 

^^^^ ^ fni(expg) _ gp -I- oi exp^ H h a„^ exp(nig) 

G„2(exp^) fco + ^iexp^H 1- exp(n20 ' 

and 

_ F„l(exp(^^)) _ qq + ai exp(ig) H h a„i exp(mi^) 

G„2 (exp(«<^)) 6o -I- bi exp(i^) H h a„2 exp(m20 ' 

where ao,--- ,a„j and 6or'' i constants to be determined, = — 1. 

Instituting the assumed solutions (4) and (5) respectively into the equation (3) 
and according to balance principle we will derive a relation of ni and n2, and 
furthermore clearing the denominators we will get two polynomials of exp ^ and 
exp(i^) respectively. In order to determine the values of parameters we set 
the coefficients of those two polynomials to be zeros, thus we get systems of 
algebraic equations. Solving the algebraic equations systems we will determine 
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those parameters and obtain the corresponding exact solutions. Of course we 
can use software such as Mathematica or maple to solve corresponding algebraic 
equations system. But our results in present paper are computed by hand. If 
computation is more complicated we will need Maple or Mathematica. 

It is obvious that our method includes the tanh-method and its generaliza- 
tions as special cases. 



3 Applications to DDFs 

Example 1. Langmiuir chains equation 

— = W„(Un+l - Wn-l), (6j 

which arise in the study of langmiuir oscillations in plasmas, population dynam- 
ics, quantum field theory and polymer science([TTl [TH [I3|). It is also named 
Volterra lattice. Under the wave transformation Un(t) = u(^), ^ = dn + ct + d., 
the Eq.(6) becomes 

cu'{0^u{OH^ + d)-u{^-d)]. (7) 

Instituting the expression (4) and (5) respectively into the Eq.(6) and according 
to the balance principle, we have rii = n2. Through detail computation, we find 
all solution are trivial solution u = constant in cases < ni = n2 < 3. In order 
to give nontrivial solutions, we take ni = n2 = 4. Thus for convenience we can 
take the solution assumed as following form 

^ A{ao + ai cxpg + a2 exp(2^) + Q3 exp(30 + exp(4^)) 

~ bo + bi cxp ^ + 62 cxp(20 + ^3 cxp(30 + exp(40 ' ^ ' 

and 

_ A{aQ + ai exp(zg) + 02 exp(2i^) + 03 exp(3zg) + exp(4ig)) 
' ~ bo + bi exp(i$) -I- 62 exp(2iO + ^3 exp(3iC) + exp(4i^) ' ^ ' 

where A,ai and bi are constants to be determined, i = 0,1,2,3. According 
to the algorithm in section 2, in real exponent function case, we obtain the 
following results: 

2c 

A = — , „ , flo = &o 7^ 0, ai = 61 = as = 63 = 0, 

exp 2d — exp(— 2d) 



a2 = (exp(2d) + exp(-2d)-l)62, ^ ^ / «o(exp(4d) + exp( 4d) 2) ^ 

y exp(2d) + exp(— 2d) — 2 

in complex exponent function case, wc have 

^ = ^T^TTT' ao = 5o 7^ 0, 01 = 61 = 03 = 63 = 0, 

sm(2d) 
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02 = (2cos(2d) - 1)62, 62 = ±2cos(d)7^. (11) 

Thus we obtain four kinds of traveling wave solutions to Langmiuir chains equa- 
tion 

_ ao + 2cosh(d)(2cosh(2d) - l)^cxp{20 + exp(40 

un[t)- csecn^dx ao + 2 cosh(d)VH^ exp(20 + exp(40 ' 

ao - 2cosh(d)(2cosh(2d) - l)>^exp(2g) + cxp(4e) 
u„(t)- csechZdx ao - 2 cosh(d)yi^exp(2e) + exp(40 ^ 



c oq ± 2 cos((i)(2 cos(2(i) — l)y^exp(2i^) + exp(4i^) 
" sin(2(i) ao ± 2 cos((i)yao exp(2i^) + exp(4i^) 

- c ,cos(20±cos(d)(2cos(2d)-l)^^^^^^^. ^^^^ 



sin(2(i) cos(2^) ± cos{d) 

and 

c ao ± 2 cos((i)(2 cos(2(i) — l)y^exp(2i^) + exp(4i^) 
" sin(2(i) ao ± 2 cos((i)yao exp(2i^) + exp(4i^) 

- c ,Bin(2e)±eos(d)(2cos(2d)-l)^^^^^^^^ ^^^^ 



sin(2d) sin(2^) ± cos(d) 

where ^ = dn + ct + S, and d, c, (5 are arbitrary constants. In special, the param- 
eter S in the Eqs.(14) and (15) have been rescaled. 

Remark 1: In the case ni — 772 > 4, the computation is more complicated, 
we will consider it in the future. 



Example 2. A general lattice equation 

ul^it) = {a + (3un{t) + 7<(i))(un+i(i) - u«-i(i)) (16) 

This general lattice equation includes some lattice equations such as Hybrid 
lattice ([2]), discretized mKdV lattice([Tll H]) as special cases. For examples, if 
a = — 1, we get Hybrid lattice 

<(i) = (1 - fiun{t) - 7<(i))K-i W - un+i{t))- (17) 

if /? = 0,7 = -1, we get discretized mKdV latticefp^ IB]) 

<W = - <W)K+i W - wn-i(<)); (18) 

Thus we consider the general lattice (16). Under the wave transformation 
Un{t) — u(^), — dn -\- ct + 5^ this equation becomes 

-u'iO = (- + -u{^) + u^imniS, +d)- d)}. (19) 

7 7 7 
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For simplicity, sometimes we make a transformation for u as follows 

/3 



then the Eq.(14) becomes 



^^'(0 = - + v'imvi^ +d)- v{i - d)}. (21) 

Instituting the assumed solution (4) and (5) into the Eq.(19) and according 
to the balance principle, we get rii < n2- If we take n2 = 1, we find there 
is only trivial solution u{^) = constant, therefore we take n2 > 2. Here we 
only consider the case n2 = 2 for simplicity, other cases will be dealt with in 
the future. When ni = 0, we get only trivial solution u = constant. When 
ni = l,for convenience we take the assumed solution as follows 

u(n - -4(ao + expO 

"^^''"6o + 6iexpC + exp(20' ^ ' 

and 

u(n = ^("o + ^^P('^)) (23) 
bo + h cMiO + cxp(22e) ' ^ ' 

where A ^ and ao,5o:^i J^re constants to be determined. According to the 
algorithm in section 2, if a = 0, then no solution is obtained, so we let a ^ 0. 
In real exponent function case we have 

c = a{expd — exp(— d)}, ao = 0, bi = 



ajexp d + exp(— c?) — 2} ' 

42 r 2! \-^\ 

, ^ ^ ta{exprf+exp(-d)-2} ^ If ^ , . 

° a{exp(2d) + exp(-2d) - 2} ' ^ ' 
in complex exponent function case we have 

c=2asm(rf), oq = 0, bi = — — 7^,bo = -p. 77-7. — rr— , (25) 

a(2cosa— 2) Q!(2 cos(2d) — 2) 

and here a and A arc arbitrary nonzero constants. Thus the corresponding 
solutions of the Eq.(16) are as follows: 

Un{t) = 

A exp(rfn + 2a sinh(rf)f + 5) ^^g^ 

^^^rfeg^'^ + ^'"""^ti^^''^'' + ^M2{dn + 2asinHd)t + S)) ' 

Unit) = 
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A cxp{i{dn + 2a sm{d)t + S)) 



^'^ 'T^r^ + + eMHdn + 2asin(d)i + S)) 

(27) 



I ^ + 



and 



n{t) = , . (28) 



-LOW 4Qain^(d/2)+^ ■ f , fl) 

=^^V 4asin^(d) Sl"^? + 4a sin^(d/2) 

Correspondingly if a 7^ 0, /? = 0, 7 = —1, we obtain the exact solutions to 
discretized mKdV lattice equation (18) as follows: 

Un{t) = ±sinh(d)Vasecli(dn + 2asmh{d)t + 6), (a > 0); (29) 

Unit) = ±sin(d)v^sec(dn + 2asin(d)t + (5), (a < 0); (30) 

Un{t) = ± sinli((i)v'— Qicsch [dn + 2asinli((i)t + S), (a < 0); (31) 

and 

Un{t) = ±sin(d)-yacsc(c!n + 2asin((i)i + (5), (a < 0), (32) 

where 5 have been rescaled. If a = —1, we obtain the exact solution to Hybrid 
lattice equation (17) as follows: 

A exp(dn — 2 sinh(d)f + 5) 



^^'S^"^ - ^"'^^tC^^W^S'^'^'^ + ^-P(2(^" - 2«inh(d)^ + <5)) 

(33) 

u^(t) = , ^ _ ^ , (34) 

Y 'l"sin^'(d) COS(0 - 4sin2^(d/2) 



and 

(i) = , ^ . (35) 

_l_9 . / ^~4sin^(d/2) • t /3 

To our knowledge, the solutions (29-35) all are new. 

When rii = 2, for convenience we take the assumed solution to the Eq.(14) 
are as follows 

_ MflQ + ai exp^ + exp(2C)) 
"(^^- 6o + &iexpC + exp(20 ' ^ ^ 

and 

^ ^(qq + «icxp(fO +cxp(2/:0) 

6o + 6iexp(iO+exp(2iO ' ^ ^ 
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where A ^ and ao, ai, bo, bi are constants to be determined. According to the 
algorithm in section 2, we have the following three families solutions: 

Family 1: for real exponent function, 

4ja - 0^ 

c = tanh(rf), fli = 6i = 0, 

27 



A = -l± tanh(.) X ^/y^, ao = ^ ; (38) 



for complex exponent function, 



47a - 

c = tan(rf), ai = 61 = 0, 

27 



/3 - 4a7 V/?^-4a7)6o 

^= -^±tan(ci) X — , 00 = — y-„ ^ (39) 

^7 ^7 /3 =F V - 4a7 

where — 4a7 > 0. Thus the corresponding solutions are as follows: 

/3 , JW^^^ -60 + exp(dn + ^i^^ tanh(rf)i + ,5) 

u„(t) = ± tanh(rf) ir^-iT2 , (40) 

^' 27 27 ^ ^ 60 + exp(dn + ^^^i^ tanh(d)t + ,5) ^ ' 

and 

/3 , - 4^7 , -60 + exp(i(rfn + ^^rf tan(d)i + ^)) 

= — tan(a) ; . (41) 

27 27 ^ ' bo + exp{i{dn + tan(d)t + 6)) 

According to the cases 60 > and 60 < 0, there are four solutions. 



Unit) = ^ ± ^^"27 "^"^ tanh(d) tanh{dn + ^"^^^ tanh(d)i + S); (42) 



Un {t) = ^± ^^"27 *anh(d) coth(dn + "^"^^^ tanh(d)t + 5) ; (43) 



u„(t) = ^ ± ^^^27 ^^''^'^^ + ^"'^^^ tan(d)t + 5); (44) 



and 



u„(t) = ^ ± ^^"27 "^"^ ^ ^^^2^ tan{d)t + 6). (45) 

From above results, if we take a 7^ 0, /? = 0, 7 = —1, then we obtain two 
exact solutions to discretized mKdV lattice equation (13) as follows: 

Un{t) = ±Vatanh(d) tanh(dn + 2atanh(d)t + S); (46) 
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Un{t) = iv^a coth(d) coth(dn + 2atanh{d)t + S). (47) 

Un{t) = ±^/atan((i)tan(dn + 2atan(d)t + (5); (48) 

and 

Un{t) — ±^/acot{d) cot{dn + 2atan(d)t + S). (49) 

If we take a = — 1, we obtain the exact solutions to Hybrid lattice equation 
(12) as follows: 



ujt) = ^ ± ^ tanh(d) tanh((in tanh( d)t + S); (50) 

27 27 27 



= ^ ± ^^^^ tanh(d) coth(dn - ^^^^^ tanh(d)t + (5); (51) 



J R2 _|_ 4^ 4-^ -|_ /32 

u„(t) = ^ ± tan(d) tan(dn tan(d)t + J); (52) 

and 

u„(t) = ^ ± ^^^^ '^'^ tan(d) cot(dn - im^{d)t + 5); (53) 

To our knowledge, the solutions (47-49) and (51-53) are new. solutions (46) and 
(50) have been given in Ref.([S]). 

Family 2: for real exponent function, 



4a7-/3^ . wn . (3 u u n 2 4(4a7-/32) . , 3, 

c= sinh(d),y4 = - — ,ao = bo,bi = 0, = — smh (d)bo; 

27 27 jj'' 

(54) 

and for complex exponent function. 



c= sin(c(),^ = -TT'^-o = Oo,0i = 0, Oi = ^ (")°o! 

27 27 /J^ 

(55) 

Hence, if 4a7 — > 0, the corresponding solutions to the Eq.(16) are as follows: 

^ (t) - ^ a± ^^'7~^^ ^i^h(d) exp(dn + sinh(d)^ + 6) 

^ 27 ^ exp(2(dn + smh{d)t + S) + 1 

= _ A(i ± ^/^"''^ ~ sinh(d) csch(dn + "^"^ ~ sinh(rf)i + 6)); (56) 
27 /3 27 



and 

V4"7-/3^ sin(d) exp(z(dn + sm{d)t + ,5)) 

27 ' exp(2i (dn + ^^^^ sin(d)t + (5) + 1 



/? \/4a7 - /32 . 4a7 - /J^ . 

= (1 ± — sin(d) csc(an H smidjt + oj); (57) 

27 p 27 
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if 4a7 — > 0, the corresponding solutions to the Eq.(16) are as follows: 

2V/3^-4a7 gj^j^^^^ g^p^^^ ^ 402^ sinh(d)t + 5) 
Unit) = —^(1 ± 1 35 ) 

27 ^ exp(2((in + ^^^^ sinh((i)t + ,5) - 1 

= _ A(i ± \/4a7-/?' sech(dn + sinh(d)t + <5)); (58) 

27 p 27 

and 



^in(i) = (1 ± S9 ■ ) 

27' e^p{2i{dn + sm{d)t + S) - 1 

= _ A(i ± \/4a7-/3^ gj^^^^ gg^^^^ ^ sin(ci)f + 5)). (59) 

27 p 27 

When a = — 1, we have the corresponding solution to hybrid lattice equation 
as follows: 

u„(t) = -ii(l ± V sinh(rf) sech(rfn + ^—t- smh{d)t + 5)); (60) 

27 p 27 

and 



Unit) = -^(1 ± ^ ^3 sin(rf) sec(dn + "^3. sin(d)i + ^)). (61) 

27 p 27 

To our knowledge the solutions (60) and (61) are new. 
Family 3: for real exponent function, 

/3(i sech(d) - 1) ± Jp'^il - \ sech(d)) - 4a7(l - sech(d)) 

^ = 7. , ao = &o, 

27 

ai = 0, 61 = -^^2^^-^^ sech(d) sinh^(d)ao, c = 2(a + + 7^^) sinh(d); (62) 
for complex exponent function, 



/?(i sec((i) - 1) ± J/32(l - i sec(d)) - 4a7(l - sec(d)) 

^ = 7, . «o = &o, 

27 

oi = 0, 6^ = - ^""^ ^^"^ sec(d) sin2(d)ao, c = 2(a + + 7^^) sin(d). (63) 



Hence we have solutions to the Eq.(16) as follows: 

Unit) = A(ao + exp(20) . ^^^^ 

ao ± ^-|^i||3sech(rf)ao sinh(d) exp(0 + exp(20 ' 
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and 

, A(ao + exp(2iC)) 

Un{t) = , ■ (65) 

'^o ± y- ^A^lfjA sec{d)ao sm{d) exp(z^) + exp(2i^) 

If ao > 0, the above solutions become 

„ A(l + exp(20) 

Un{t) = . ; (66) 

1 ± ^-^1^ sech(d) sinh(d) exp(0 + exp(2e) 

and 



1± 



-i^ig^ sec(rf) sin(rf) exp(iO + exp(2iO 



If ao < 0, the above solutions become 

u.it) = ^(-l+exp(2e)) , ^^3^ 

-1 ± i/^l^ll^ sech(d) sinh(d) exp(0 + exp(20 

and 

u.it) = , A(-l + exp(2^0) ^ ^^^^ 

-1 ± ^ sec(rf) sin(rf) exp(iO + exp(2iO 

Furthermore, when a^O, /J = 0, 7 = — 1, respectively we have the solutions 
to the discretized KdV equation as follows: 

u^it) ^ , , -4(l + cxp(20) ^^^^ 

1± ^y/-a sech(d) sinh((i) exp(^) + exp(2^) 

where asech(d) < 0. And 

Unit) = , , M-l + eM20) 

— 1 ± ^ V a sech(d) sinh((/) exp(^) + exp(2^) 



where asech((i) > 0. 



Unit) = \ , (72) 

cos^± ^y— asec(a) sm(a) 

where asec(d) < 0. And 

Asin£ 

sm^ ± ■J-^/— a sec (a) sm(a) 

where a sec(d) < 0. 

When a = —1, we have the solutions to hybrid Lattice equation as follows: 

A(l + cxp(20) 

Un{t) = . ; (74) 

1 ± sech((i) sinh((i) exp(0 + exp(20 



10 



Un{t) = , ^ = ; (75) 

-1 ± ^J^^^ sech{d) sinh(d) exp(0 + exp(20 

r 2 A COS f 
Unit) = , ; (76) 

2cosC ± ^J-^^^ sec(rf) sin(rf) 

Un{t) = , . (77) 

2sinC ± sj-^^^ sec(rf) sin(rf) 

To our knowledge all these solutions given in family 3 are new. 

4 Discussions and conclusions 

We proposed a new method named rational expansion method of exponent 
function to solve exact solution of nonlinear differential-difference equations. 
This method generalized Tanh-method and some other similar methods. Using 
our method, abundant exact solutions to Langmiuir lattice equation, Hybrid 
lattice equation and discrcmizcd mKdV lattice equation. Among them, many 
solutions arc new. Because of using the property of exponent function 

exp(^ + d) = exp ^ exp d, exp(i(^ + d)) = exp(i^) exp(id), (78) 

we transfer the addition into the multiplication, so corresponding computation 
becomes relative simple. For a lot of DDEs such as Toda lattice, discrete KdV 
lattice and Ablowitz-Ladik lattice and so on, our method is also applicable. Of 
course our method can be applied to lattice equations system. 

We must point out that there are also some problems for our method to 
study further. An open problem is listed in the following: 

Open problem: whether or not our method can always give new solutions 
when the orders of polynomials F and G in the Eqs.(4) and (5) increase step 
by step. 

Here I have an example to illustrate above problem. We consider the rela- 
tivistic Toda lattice system 

Wn(*) = (1 + OlUn){Vn - ^n-l), (79) 
'"'n{t) = Vn{Un+l - W„ + UVn+l - aVn-l)- (80) 

Take a transformation = — — we have 

u'nit) = {Un^-){Un-l- Un). (81) 

a 

Under the traveling wave transformation Un{t) = u{£,), ^ = dn + ct + S, we have 

cu'{0 = HO + -)iu{^-d)-u{0). (82) 
a 
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Instituting the Eq.(4) into the Eg. ([5^ and according to the balance principle 
yield ni =712- Moreover we can easily prove the following conclusion: 
For any positive integer fc, the equation (81) has solutions as follows: 

_ A{ao + exp(fcd + fcc + kS)) 
''"^^^ ba + cMkd + kc+kS) ' ^ ' 

where A = ^^p/ld)-! - h^-o ^ {(A+-)exp(fcd)--}fcn . ^ ^ arbitrary con- 
stants. For the complex case, we have a similar result. 

From above conclusion, although k is arbitrary, if we rescale d, c and 6, it 
is easy to see that we can't give new solutions. This result implicates that our 
open problem is difficult to answer. 
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